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O ■ ON PREDICTING THE POLARIZATION OF LOW 

O FREQUENCY EMISSION BY DIFFUSE INTERSTELLAR DUST 



P. G. Martin^ 



Abstract. Several of the current and next-generation cosmic microwave 
background (CMB) experiments have polarimetric capabihty, promis- 
ing to add to the finesse of precision cosmology. One of the contami- 
nating Galactic foregrounds is thermal emission by dust. Since optical 
interstellar polarization is commonly seen, from differential extinction 
by aligned aspherical dust particles, it is expected that thermal emis- 
sion from these grains will be polarized. Indeed, in the Galactic plane 
and in dark (molecular) clouds, dust emission in the infrared and sub- 
millimetre has been measured to be polarized. It seems likely that the 
faint diffuse cirrus emission, of more relevance to CMB experiments, 
will be polarized too. We discuss how well the amount of polariza- 
tion of this component can be predicted, making use of what is known 
about optical (and infrared and ultraviolet) interstellar polarization 
and extinction. Some constraints on the alignment of the carrier of the 
dust-correlated anomalous microwave emission can be made as well. 
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1 Introduction: CMB Polarization and Cirrus Contamination 

• I— I 

■ Several of the current and next-generation cosmic microwave background (CMB) 
] experiments have the capability of measuring linear polarization, a complemen- 
. - - 1 tary probe for precision cosmology (e.g., Laureijs 20CJH^; Polenta et aLUD05, Benoit 
et al. 121)1)41 Prunet l!!UU6|l . Among these experiments are Archeops (Benoit et al. 
IMTljl . BOOMERanG (B2K in 2003; Polenta et al. Netterfiled 2006 ) , and 

the Planck Surveyor (Laureijs "2006! Prunet '2006' Lammare et al. I2003fl . All-sky 
surveys by IRAS from 12 to 100 /xm revealed faint diffuse emission ("cirrus") ev- 
erywhere, even at high Galactic latitudes (Low et aZ. [1984). One of the diffuse 
foregrounds contaminating the CMB signal is thermal emission by this diffuse in- 
terstellar dust (e.g., de Oliveira-Costa et al. 120041 Tucci et al. I2005p . expected to 
correspond most closely spatially to the 100 /Ltm emission (S I2.1|l . This becomes 
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Fig. 1. Left: Foreground components contaminating tlie CMB at two representative 
latitudes. From www.planck.fr/headingl36.html Giard & Lagache. Right: Alternative 
presentation from de Oliveira-Costa et al. I2UU4II showing the anomalous emission as well. 

the dominant foreground above about 100 GHz (A < 3 mm = 3000 /im), and as 
the CMB fades toward higher frequencies, the rising dust spectrum becomes the 
dominant signal (Fig. QJ. Where polarization is concerned, some further caution 
is needed in defining where the dust dominates. The frequency above which the 
polarized dust intensity exceeds the polarized synchrotron background is probably 
somewhat higher because the synchrotron component is likely to have a higher 
degree of polarization. Of course, it is the spatial fluctuations, which will vary 
from component to component, that are in the end important to CMB contami- 
nation. Multifrequency observations through the transitional range are therefore 
a prerequisite for separating these signals (component separation, not a specific 
topic of this paper; see e.g., Jaffe et al. I2()()4|l . 

Some CMB experiments like B2K minimize cirrus contamination by concen- 
trating measurements in regions of low foreground column density (Netterfield 
I2006|l . One strategy for cirrus mitigation could therefore be to mask out re- 
gions of bright cirrus. But only 20% of the sky has an HI column density below 
10^° cm— 2 (Miville-Deschenes, private communication) and even that produces 
a non- negligible foreground 1 MJy sr— 1 at 100 /Ltm). The precision sought 
by instruments like the Planck Surveyor depends on broad sky coverage (Prunet 
1200611 ■ and of course Planck will in any case produce data for the entire sky. Thus, 
as an alternative strategy, one needs to measure the properties of cirrus at higher 
frequencies where the CMB is not important, and then extrapolate to lower fre- 
quencies where one does have to address component separation. To inform/confirm 
the extrapolation, the cirrus properties at the lower frequencies can be assessed 
empirically in regions of brighter cirrus. The spectral properties of cirrus are of 
course an important diagnostic of the interstellar dust. 

Since optical polarization is commonly seen, from differential extinction by 
aligned aspherical dust particles, it is expected that thermal emission from these 
grains will be polarized (Stein IT!)67(I . Indeed, in the Galactic plane and in dark 
(molecular) clouds, dust emission in the infrared and submillimetre has been mea- 
sured to be polarized fVaillancourt .2006. Hildebrand et a/. .2000i Vaillancourt et al. 
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11003 • In these dense regions, the interpretation of the polarization is complicated 
by many issues: beam dilution, distortions in the magnetic field topology, changes 
in the degree of alignment within clouds, grain evolution, and a range of grain 
temperatures and optical depth which affect which grains dominate the emission 
and polarization in various parts of the infrared to submillimetre spectrum. These 
are not considered here. 

It seems likely that the faint diffuse cirrus emission will be polarized too, and 
this will not be affected by the above-mentioned complications. Diffuse dust po- 
larization is of more relevance to CMB experiments. Polarimetric observations 
in the transitional high frequency range (see Table ^ were made by Archeops at 
353 GHz (Benoit et al. EUMjl and B2K and MAXIPOL at lower frequencies (Po- 
lenta et al. 120051 Johnson et al. I2003|) . precursors to Planck HFI. Archeops has 
reported detecting the dust polarization (§ refsobs). From the point of view of 
best examining the dust polarization, Planck has no polarimetric capability at 545 
and 857 GHz where the dust polarized intensity would be the strongest. However, 
balloon-borne experiments like PILOT (Bernard I2005|l and BLAST-pol will map 
the diffuse polarization with great sensitivity at higher frequencies. 

This paper attempts to quantify the degree to which the diffuse dust emission 
is polarized, a contribution toward both component separation and dust physics. 

WMAP (Kogut et al. I2003 >l and Planck have polarimetric capability at lower 
frequencies, down to 23 GHz and 30 GHz respectively, which can be used to 
examine the polarization of two other CMB foregrounds, synchrotron emission and 
the dust-correlated anomalous microwave (hereafter simply "anomalous" ) emission 
(Draine & Lazarian J^998 : Finkbeiner et al. 119991 Lazarian & Finkbeiner 20031 
Finkbeiner IWil Finkbeiner et al. IMTH de Oliveira-Costa et al. EITOI Davies 
I2()()()l Davis 1^)06(1 . Some constraints on the alignment of a dust carrier of the 
latter emission, which would be diagnostic of the emission mechanism, are made 



Table 1. Angular Resolution^ of High- Frequency CMB Polarization Experiments 



V (GHz) 100 143 150 217 240 345 353 545 1250 

A (^m) 3300 2100 2000 1380 1250 870 850 550 240 



Archeops 13 
BOOMERanG 2K 10 7 7 

MAXIPOL 10 10 

PILOT 2.3 1 

Planck HFI 9.2 7.1 5 5 



^ FWHM in arc minutes 
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here as well 15. 5|) . 

1.1 A Legacy from Optical Polarization 

Both optical interstellar polarization and polarization of thermal emission depend 
on aspherical (flattened or elongated) grains that are aligned. The orientation of 
the electric vector of the emitted radiation is along the long axis of the mean grain 
profile projected on the plane of the sky, while that of the optical polarization, 
being of transmitted light suffering differential extinction, is orthogonal to this. 
Particular ingredients in alignment theories are that grains are rapidly spinning 
about their axis of largest moment of inertia and that this is on average aligned 
with the Galactic magnetic field, so that consistent with observations the electric 
vector of optical polarization is along the direction of the projected magnetic field, 
while it is perpendicular for emission (see Fig. EJ. There are several potentially 
viable alignment mechanisms as discussed elsewhere (Lazarian & Finkbeiner 120081 
Lazarian & Yan '2004 Roberge 2004'; Vaillancourt 200B|) and in the Appendix. 
Nevertheless, it is still difficult to predict the degree of alignment with certainty 
ab initio, or as will become clear, gauge other important factors, and so a focus 
of this paper is to assess how well the degree of polarization of the diffuse cirrus 
component can be predicted semi-empirically. 

To predict the degree of polarization of the dust-related CMB foreground(s), 
we draw on what can be deduced about alignment from optical (and infrared 
and ultraviolet) interstellar polarization and extinction, concentrating on diffuse 
dust associated with atomic gas. Still, it should be cautioned that the optical 
measurements are by necessity largely for column densities with Eb-v > 0.1 
and so Nh > 6 x 10^° cm~^. This corresponds to lines of sight with 100-/j,m 
cirrus brightness /loo > 6 MJy sr^^. While there will be many such directions 
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Fig. 2. The relationship of the polarization of emitted and transmitted radiation to 
the mean grain profile, aligned with respect to the Galactic magnetic field. In the 
far- infrared and submillimetre, there is more radiation emitted with polarization E- 
vector along the long axis. This is also the axis along which there is more extinction, 
and so the E-vector of the net transmitted radiation is orthogonal to this axis, and to 
the emitted polarization and the magnetic field. Graphic from Ponthieu & Lagache at 
www.planck.fr/article263.html. 
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in which to examine the polarization of the low-frequency foreground cirrus with 
forthcoming experiments, for CMB applications one needs to work at even lower 
column densities too. Another issue for any intercomparison of optical polarization 
and the polarization of low-frequency dust emission is that the Galaxy is optically 
thin for the latter. Therefore, we see the whole Galaxy, or right out of it, unlike 
probes with stars which rely on differential extinction along that path. But at high 
latitude, most relevant to the CMB, the effective paths might not be too dissimilar 
if the stars used were sufficiently distant. 

A lot of attention goes into the degree of polarization p — P/I, the ratio 
of the polarized intensity P to the total emission /. Non-polarizing grains (or 
other emission components entirely) can contribute to /, potentially confusing the 
interpretation of any spectral dependence of p. We therefore also emphasize the 
importance of examining P and its spectral dependence directly, as a diagnostic of 
the aligned grains. P is of course what is relevant to the calculation of the power 
spectrum of the E and B modes of the polarized intensity. We do not attempt to 
predict these (Prunet et al. 119981 Tucci et al. I2005|l since they depend in addition 
on the spatial variation of the degree and orientation of the alignment, whose 
statistical properties are not known (i.e., the statistics of P are different than the 
statistics of Pem-^)- However, the dust contribution to the power spectrum should 
scale with the same frequency dependence as P^ (from dust), unless hoped-for 
"spatial - electromagnetic frequency decoupling" breaks down, i.e., there were 
different populations of grains with different alignment statistics whose relative 
contributions to P (and /) changed with frequency. A TE measurement of the 
dust emission at 353 GHz by Archeops (Ponthieu et al. I2()()5|l provides a basis 
from which the potential contamination at lower frequencies can be estimated 
more directly. 

2 Polarized Emission 

The polarization of thermal emission by dust was first noted and estimated by 
Stein H1967|l . Hildebrand l)1988|l has since made more detailed calculations re- 
garding infrared polarization and Hildebrand & Dragovan l|1995|l have analyzed 
far-infrared observations. We begin our estimate on the same basis, from electro- 
magnetic scattering calculations of the cross-section for emission (equal to that for 
absorption), Cem, for single grains (for much of the rest of this section we suppress 
the subscript em). What is often computed for a single grain (e.g., Martin . 1975|I 
is an efficiency factor Q — C /A, where A is some measure of the area of the grain 
(whether projected or actual does not matter here). (Note that this is not Stokes 
Q.) But first we recall one of the challenges, that there are many different grain 
components potentially involved. 

2.1 Grain Components: Lessons from IRAS and ISO 

Interstellar extinction suggests a range of grain sizes and compositions f S I3.1|l and 
the infrared spectrum of cirrus in the range 3 to 1000 /Ltm is interpreted as the sum 
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Fig. 3. Spectral components of high Galactic latitude cirrus in the infrared and im- 
plications for grain sizes. Left: Desert et al. lflW))l . Right: Li & Draine PM) . See 
text. 



of several components too (Fig.O. 

Near 100 /j,m and longer, the range of interest here, the emission is thermal 
emission by "big" grains (size ~ 0.1 /im) in equilibrium with the local interstellar 
radiation field. 

The mid-infrared emission at 60 and 25 /im is in excess of what can be at- 
tributed to the big grains, and is interpreted as non-equilibrium emission by grains 
that are very much smaller (VSGs: 0.007 /im = 70 A = 7 nm), such that a single 
absorbed ultraviolet photon raises the temperature momentarily above the equi- 
librium value. 

The 12 /im and shorter emission has many spectral features indicative of poly- 
cyclic aromatic hydrocarbons (PAHs). This can be thought of as either non- 
equilibrium emission by tiny grains/PAHs (1 nm) or the result of internal conver- 
sion in large molecules (see e.g., Draine & Li l2UUlfl . 

All of these smaller grain components of course radiate at longer wavelengths 
too. But because the big grains contain most of the mass, it is they that dominate 
in the submillimetre. 

Tiny grains (including PAHs) also spin rapidly and emit microwave radiation 
which could be another foreground contaminant of the polarized CMB, if aligned. 



2.2 Submillimetre Spectrum: Total Intensity from an Ensemble of Grains 

The emitted total intensity / summed over grain populations with column density 
Ng and temperature Tg, and acknowledging all "other" components (CMB, free- 
free, synchrotron, anomalous emission, point sources, etc.) is 

C (T) 

h = B,{T)NAQ,{T) = B,{T)NV^^ + hther- (2.1) 

Here C is the mean of the two effective cross sections for orientations of the electric 
vector parallel and perpendicular to the long axis of the projected mean profile of 
the aligned grains ("par" and "per" below in eauation l2.2|l . Normally C depends 
on the grain size (characterized by some dimension a) but in the millimetre to 
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submillimetre range of interest, a is much smaller than the wavelength A and 
the absorption cross sections per unit volume V, C/V = Q/a, is independent of 
a. The contributions to I (and P) are therefore weighted by the volume; in the 
diffuse interstellar medium, large grains carry most of the volume (or mass) . It is 
often a good approximation in this limit that C/V = Q/a = av^ , where fi is the 
spectral index of the emissivity (e.g., (3^2 for many materials) and a is the scaled 
opacity (with the basic frequency dependence removed), that varies from material 
to material. The thermal emission spectrum of each separate grain component is 
then a modified greybody emission, ~ By{T)y^ , with the appropriate T and (3. 

If there are multiple grain components with different properties (including 
alignment), then the contributions to / (and P) are weighted by the volume, 
a, and the Planck function at the appropriate temperatures, potentially making 
/ deviate from a simple modified greybody with a frequency- independent /?; nor 
would (3 be necessarily representative (diagnostic) of any of the grain materials. 
Spectral changes would be most pronounced near the peak of any of the Planck 
functions, because of changing relative weights. 

The wider the frequency range, the less likely it is that the opacity can be 
described by a single (3. Changes in (3 with frequency over the range of interest 
are seen in laboratory experiments on amorphous silicates (Boudet et al. I2()()5|l . 
In fitting the spectrum, this provides an alternative to invoking an extra diffuse 
cold (< 10 K) dust component (Finkbeiner et al. 119991 Bourdin et al. l2UU2f) : 
although this provides a good fit and interpolation formula, one must be wary 
about taking the model literally, since dust that cold in the local diffuse ISM seems 
unphysical. Also there is the possibility, often ignored, that both a and (3 (thus the 
opacity) could be temperature dependent, again supported by the same laboratory 
measurements on amorphous materials, where the effects are attributed to intrinsic 
processes (two level systems, disordered charge distribution) in a single amorphous 
(silicate) material. This might also help explain why (3 appears to change with T 
(Boudet et aL 12(1051 Dupac et a,l. \2()()[M . though here it can be noted that this is for 
T higher than in regions of relevance to diffuse CMB foregrounds. As discussed in 
§ 15.41 an important implication might be that silicates provide a more dominant 
contribution to the emission at lower frequencies (§ 15. 1|) , of relevance to polarized 
emission (and its detailed frequency dependence) since, among grain components, 
large silicates at least are aligned I3.2|l . 

2.3 Polarized Intensity 

The emitted polarized intensity P (Stokes Q in the appropriately-rotated coor- 
dinate system) summed over all (aligned) grains, and including a shorter list of 
"other" polarized components, is similar: 

Peva — Piy{Pg)-^g'^{Cpar Cp^^^ + Pother- (^■^) 

Since Cpar > Cpen the emitted light is polarized with electric vector parallel to 
the long axis of the mean grain profile projected on the plane of the sky. The 
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minus sign invokes a convention, a reminder of the orthogonal orientation of the 
electric vector for emission relative to that of the transmitted light (see the sign 
convention in equation 13. 

2.4 Emission Cross Sections for "Small" Particles 

In the relevant limit in which a is much smaller than the wavelength A simple 
analytical formulae can often be used. These are particularly straightforward to 
calculate for spheroids (or ellipsoids for that matter) that are homogeneous. In 
this limit there are only two principal cross sections to evaluate, corresponding to 
electric vector orientation parallel and perpendicular to the symmetry axis: 

a V c ^L||,^(m2 - 1) + 1^ ^ ' 

The principal cross sections are analytic functions of the complex refractive index 
m (hence dependent on grain material) and the "polarizability" factors L which 
depend on the shape and axial ratio (e.g., Martin [l975fl . 

Interstellar grains are spinning and the largest amount of polarized emission P 
would arise for what we will call perfect spinning alignment (PSA), corresponding 
to rotation about the grain axis of largest moment of inertia, with all spin axes 
aligned along a common axis in space (the interstellar magnetic field) which in 
turn is in the plane of the sky. 

For spheroids with perfect spinning alignment (PSA, as opposed to picket fence 
alignment, PF) we have Qpar,per = Q±,\\ for oblate shapes and Qpar,per = iQ\\,± + 
Q±)/2 for prolate shapes. 

2.5 Degree of Polarization of Aligned Grains paem 

For a single grain population (size, composition, and orientation, leading to the 
same temperature Tg), the dependences on NgV and B^iTg) would cancel out in 
the ratio P/I and so 

Qpar Qper , 

paem = • (2.4) 

^par ~r ^per 

In this simple situation, the degree of polarization is a property characteristic of 
the grains (composition, shape, alignment), but not their column density. 

For m appropriate to "astronomical silicate" (Draine & Lee 119841 Rouleau & 
Martin ri991|l . Figure ^ illustrates paem calculated for several axial ratios. Also 
shown is one curve for an amorphous carbon material (ACl; Rouleau & Martin 
[TMT|i . 

There are several things to note. The degree of polarization is potentially very 
large, even from particles that are only mildly aspherical (axial ratio < 2). The 
value of puem depends on the shape, axial ratio, and degree of alignment, and so 
by itself is not diagnostic of the specific grain material. Even though Q/a has a 
strong spectral dependence as described above, this too cancels out in pa^^rn- There 
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Fig. 4. pflem for perfect spinning alignment, as a function of shape and composition. 
Silicate: oblate with axial ratio 1.4 (dash) and 2 (solid); prolate with axial ratio 2 (dash- 
dot). Amorphous carbon: oblate with axial ratio 2 (dots). 

can be a slight frequency dependence of pOem across the range of interest, due to 
changes in m. 

In contrast, the spectral dependence of the observable P does retain the i^^ 
dependence. The wavelength dependence of optical polarization indicates that 
only large grains are preferentially aligned, which reinforces their dominance in 
the volume-weighted P and makes it more likely that the spectral dependence of 
P could be a useful diagnostic of the material of the aligned grains. For example, 
if over the range ~ 100 - 1500 GHz interstellar amorphous silicates had spectral 
index variations as measured by Boudet et al. (j2U05|l or inferred by Li & Draine 
H2001|l and these were the aligned grains (silicates certainly are aligned; see 13.2(1 . 
then this spectral behaviour would be imprinted directly on P with little effect on 
pflem (see the illustration in I5.4|l . 

Non-aligned grains and lother dilute the degree of polarization that would be 
produced by the aligned grains. The spectral dependence of P would not change, 
but the interpretation of would become more difficult. 

More generally still, spectral dependence in the net Pem can be introduced by 
frequency-dependent relative weighting in P and/or /, from different degrees of 
alignment, and/or temperatures Tg and opacities (many models implicitly imply 
these effects, e.g., Finkbeiner et al.^M Bourdin et a/. 12003 Li & DraineEOOH), 
and from the contributions Pother and lother- 

2. 6 Challenges 

As mentioned above, there is a wide range of grain sizes and compositions. In order 
to predict the polarized emission, we need to know which of these grains are aligned 
and how well. This is wrapped up with grain shape too: how flattened/elongated 
are the particles? In summary, pem (and P) depends sensitively on axial ratio. 
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grain orientation, and composition, none of which arc well known. 

3 Optical (and Infrared and Ultraviolet) Interstellar Polarization 

Fortunately, though the challenges can be met; the appropriate combinations of 
these unknowns can be constrained using interstellar polarization in the optical. 
First we reinforce one of the challenges, the multiple components suggested by 
interstellar extinction curve, T^xi^^)- 

3.1 Grain Components: Lessons from the Extinction Curve 

Interstellar extinction is caused by scattering plus absorption by grains comparable 
in size to the wavelength. The dependence of the cross sections Cex on size and 
frequency cannot be obtained by the small particle approximation described above 
and require instead complex electromagnetic scattering calculations such as by 
Mie theory. The grain models are used to interpret observations of the interstellar 
extinction curve (Fig. [Sj left). 

One basic conclusion is that the continued rise in extinction into the ultraviolet 
requires smaller and smaller grains, a range already known prior to the interpre- 
tation of the infrared emission. There is a well-characterized "bump" at 2175 A 
thought to arise in carbonaceous grains. This is another indication of separate 
grain components. In the infrared at 10 /im (Fig.jSJ right), there is a prominent 
silicate absorption feature (Si-0 stretch), the strength and shape of which requires 
that most of the available Si is depleted in amorphous silicate grains. The near- 



Fig. 5. Frequency dependence of interstellar extinction. Left: Extinction curves in the 
optical and ultraviolet showing a range of behaviour as a functon of Rv (Cardelli et 
a/. ©89 ); in the diffuse medium Rv ~ 3.1. The sole spectral feature is the prominent 
"bump" at 2175 A. Right: In the infrared there is a power-law decrease, plus a distinctive 
10-^m silicate feature (Martin & Whittet HMTl . 




Computed 
Observed 



BOt56 5£4 
R-2.75 



p. G. Martin: On the Low Frequency Polarization of Cirrus 11 



infrared extinction follows a power law with index near 2 (Martin & Whittet ll990|l 
which can be modeled as the scattering by large grains comparable in size to the 
wavelength. 

In summary, grains come in many sizes (perhaps a function of composition). 
This raises many questions. Which grains produce the extinction in the optical 
and ultraviolet? Which grains produce the submillimetre emission? Which grains 
polarize in the optical and ultraviolet? Does this result in significant submillimetre 
polarization? 

3.2 Lessons from the Polarization Curve 

At least some of the grains are aspherical and aligned, so that in the plane of 
the sky the ensemble average grain profile is elongated, empirically with the long 
axis of this profile oriented perpendicular to the direction of the magnetic field 
B. Differential extinction, according to the orientation of the electric vector E 
with respect to the mean projected grain profile, produces a net polarization of 
transmitted light. Since there is greater extinction for E parallel to long axis, then 
E of the transmitted radiation is parallel to short axis, hence parallel to (Fig. [21). 

Empirically, the "polarization curve" , the frequency dependent degree of polar- 
ization Pf,x{i'), increases in the infrared to a peak in the optical and then declines 
into the ultraviolet (Fig. El left). Computations for single sizes indeed show a 
decline in polarization before the peak is reached in extinction (Roger & yi&v- 
t in 11979(1 . But quantitatively this plus the continued rise of extinction into the 
ultraviolet, which implies smaller and smaller particles, plus the low ultraviolet 




Fig. 6. Frequency dependence of interstellar polarization. Left: Polarization curve from 
the near-infrared to the ultraviolet, normalized to the maximum polarization in the vi- 
sual (Martin et al. I1999II . There is no strong ultraviolet polarization "bump" as in the 
extinction curve. Right: In the infrared there is a power-law decrease, plus a distinctive 
polarized 10- /^m silicate feature (Martin & Whittet I199UII . In lines of sight to an em- 
bedded source, such as this toward the Becklin-Neugebauer object in OMC 1, there is 
polarization at S.l^tm ice band, interpreted as a thin frost on the aligned silicates. 
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polarization, leads to the important conclusion that only the larger grains are as- 
pherical and aligned; the aligned grain mass distribution is very deficient in small 
particles (Kim & Martin ,199.5j) . 

The strong 2175 A bump does not have a correspondingly dramatic polarization 
feature at 2175 A. A weak feature is seen only in two of the 28 lines of sight in 
the Galaxy observed by WUPPE and HST (FOS), despite having sufhcient S/N 
(Martin et al. I1999|l . It is not a common phenomenon. The polarization feature is 
at the same position as the extinction bump, has a positive excursion, and shows 
no change in position angle. But it is at least two orders of magnitude smaller 
than the theoretical maximum for perfectly aligned graphite carriers (Martin et 
al. I1995| . Thus, either the alignment is quite incomplete or only a small fraction 
of the grains is aligned. The result that the (presumed small) particles giving rise 
to the 2175 A extinction bump are very poorly aligned is consistent with the poor 
alignment of the (other) small particles which give rise to the continuum ultraviolet 
extinction. A small residual alignment might be arise from the Davis-Greenstein 
process (Wolff et flZ. jTMTD . 

Significantly, the feature at 10 /zm is polarized (Fig.|Hl right), indicating that at 
least the silicate grains are aspherical and aligned. Details of the relative changes 
of Pex and Tex across the feature constrain the band strength and the shape and 
axial ratio; Hildebrand & Dragovan H1995|l conclude that the silicate particles 
are oblate with axial ratio ~ 1.5. It might be relevant that this is not unlike 
the individual sub-grains found in fluffy silicate agglomerate interplanetary dust 
particles, and the further insight from cometary material collected by Stardust is 
eagerly anticipated. 

For embedded sources in molecular clouds, the polarized ice feature at 3.1 /xm 
can be interpreted as a thin coating on the aligned silicates (as predicted: Martin 
I1975|l . The rise in the infrared is a power law (Martin & Whittet ll99UI Martin 
et al. I1992|l not unlike that for extinction, which again can be modeled as the 
(differential) scattering by large grains comparable in size to the wavelength (Kim 
& Martin [TM^ . 

In summary, only the larger grains are aspherical and aligned. The 10 /im 
polarization shows that the silicate component is aligned, and the (oblate) parti- 
cles have an axial ratio which is modest (near 1.5). It seems at least consistent 
(certainly adequate) to model the polarization with silicates alone. 



3.3 Modeling 

The observable Pex{i^) increases with the column density of aligned grains. It is 
thus useful to normalize this out by considering the ratio with respect to extinc- 
tion Tex.v which also builds up with column density. Suppressing the frequency 
dependence of the cross sections, 

/Pn _ Ngi^jCpar — C per) ex ,„ 

(see, e.g., Martin [r974|l . This has a form that is deceptively like that for pem, apart 
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Fig. 7. Typical frequency dependence of (p/T)ex in the diffuse interstellar medium (solid 
curve). This is derived by dividing a typical polarization curve (dashed; Martin et al. 
I1999II by a typical extinction curve (dash-dot; Cardelli et al. 1989 for Rv = 3.1, normal- 
ized to TV = 1, with scale on right). Both p^x and (p/T)ex are normalized (scale on left) 
according to the maximum value observed at V (see Fig.|5Jl. The ratio {p/T)ex at 2.2 ^m 
(K band, left vertical dashed line) is significantly higher. 

from the weighting of the latter with the Planck function (see equations 12.11 and 
12.21 and 12. 4|) . However, the further simplifications arising from the approxima- 
tions for Qem for small particles (e.g., simple volume weighting; size-independent 
behaviour) do not apply. 

Figure shows the very different frequency dependences of t^x (dash-dot) and 
Pex (dash) derived from observations in the diffuse interstellar medium (??ccm, 
??mar99). Both Tf,x and Pex are sensitive to the actual size distributions, for 
all the grains and the subset that are aligned, respectively. Therefore, a self- 
consistent model to predict {plT)ex and eventually Pem requires size distributions 
that reproduce the frequency dependence of both the polarization and extinction 
curves. 

4 Constraining Alignment and Shape 

Both the degree of alignment and shape influence the degree of polarization and 
so have to be constrained in order to predict Pem ■ 

4.1 Imperfect Alignment 

For the ensemble of grains, disalignment can be characterized by disorientation of 
the direction of the magnetic field along the line of sight, of the angular momentum 
vectors of the spinning grains with respect to the magnetic field, and of the spin 
axes (or body axes) with respect to the angular momentum vectors (e.g., Greenberg 
[1968. PurceU & Spitzer .1971. Martin . 1975. Hildebrand .1988|i) . These affect the 
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emission P and /, and Pex and t^x, in subtly different ways. 

For grains small compared to the wavelength, as appropriate for pem, the cross 
sections for arbitrary orientations can be written as linear combinations of the 
principal cross sections, so that the average P and I over a distribution functions 
can be obtained more straightforwardly. For P the effects are actually separable, 
and not functions of shape, size, or composition, which can be useful; Greenberg 
gives a size-independent "Rayleigh reduction factor" . Of course, in the 
general case, the appropriate sums in P over the grain components, which might 
have different alignment, still have to be performed. In the optical there are no 
such approximations, but there is some indication from representative calculations 
(e.g., ??green, ??rogers) that for the same disalignment the fractional reduction in 
Pex is quite similar to that of P. Further calculations to quantify this would be of 
interest, but are beyond the scope of this paper. 

There are related effects through the cross sections, which are not separable 
even for small particles, that tend to increase the emission / and Tex relative to 
perfectly aligned grains, further reducing pem and {p/T)ex, respectively. The ef- 
fects are small for small axial ratios, but can be more substantial as grain flattening 
or elongation increases. As discussed below, in the interstellar medium even the 
maximum optical {p/T)ex corresponds to either axial ratios close to one, or close to 
random orientation for more elongated particles, with some independent evidence 
for the former. Thus for / and Tex calculations using spherical particles and/or 
random orientation might be close to appropriate. But without complete calcula- 
tions for disaligned grains, this is one major source of systematic uncertainty in 
the applying "recipe" below 14. 

Relative to perfect spinning alignment (PSA), which for computational sim- 
plicity is often used, one can define a "reduction factor" R which quantifies how 
the polarization is reduced when the grains are less than perfectly aligned. From 
the discussion above. Rex for (j>/T)ex could be different than Rem for Pem and be 
frequency dependent. However, they are likely to be quite similar, and because 
little is actually known about the details of alignment. Rem — Rex is a pragmatic 
approximation to adopt. 

4.2 An Empirical Constraint 

The observed amount of optical polarization per unit extinction provides an empir- 
ical measure of the asphericity and degree of alignment, whence R for a particular 
grain model with specified shapes and sizes. 

FigureElshows {j)/T)ex derived from the empirical curves for Tex and Pex- There 
is a strong systematic spectral dependence. The normalization is for lines of sight 
in the interstellar medium for which the most optimal combination of shape and 
alignment is achieved (Fig. IHJ. This corresponds to a value {p/T)ex — 0.0267 at 
the V passband or about 0.06 at K. 

Hildebrand & Dragovan H1995|l adopt a maximum, p/r = 0.065 at the K band. 
While this is from lines of sight with higher extinct ion, and perhaps higher density, 
observed by Jones et al. I)1992|l . and so perhaps less appropriate to the cirrus 
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Fig. 8. Observed values of optical polarization near V plotted against colour excess Eb-v 
(Serkowski et al. 119751 . For a typical ratio of total to selective extinction Ry ~ 3.1, the 
upper envelope corresponds to {p/t)^^ = 0.0267 at V or about 0.06 at the K band 
(Fig. |7|l. For directions in which the magnetic field is more along the line of sight, and/or 
for which there is a lower degree of alignment, (p/r)ei is lower, typically by a factor two, 
and Pern of the low frequency emission will be lower too. 

application here, the value is nevertheless very similar. 
4.3 A Bootstrapping Recipe 

For their analysis of far-infrared polarization in more dense regions of the inter- 
stellar medium (clouds), complicated regions for which a maximum pem ~ 9 % 
is observed (Vaillancourt I2U(J6 ^ Hildebrand et al. I2(J(J(J|I , Hildebrand & Dragovan 
P995|l gauge the effect of disalignmcnt by first comparing pem for their model at 
2.2 /im with {p/T)ex observed at 2.2 fim. In practice, there are two problems with 
their approach. First, they assumed that the formulae involving cross-sections 
for pure absorption would apply to both pem in the far-infrared and {p/T)ex at 
K. However, the slope of the observed power-law dependence of extinction and 
polarization near K indicates that this involves grains of size comparable to the 
wavelength, where scattering is very important. The adopted formulae do not 
really apply at 2.2 /im. Second, even if this were the case, the evaluation of Pe 
near 2.2 /xmfor silicates is very sensitive to how "dirty" they are, as quantified 
by the size of k, the imaginary part of the refractive index in the near- infrared 
and optical. This is poorly known. On the other hand, how "dirty" the silicates 
are is not relevant to the low frequency polarization. Therefore, for silicates, it is 
hard to bootstrap from the near-infrared to the far-infrared or submillimetre with 
confidence. Addressing the first problem with more realistic scattering models also 
eliminates the second. 

To proceed, as stated in S 13.31 we must compute (p/T)ex according to equa- 
tion using realistic size distributions and an electromagnetic scattering theory 
(Mie-like) for particles comparable in size to the wavelength. This is not impos- 
sible, and so a practical bootstrapping recipe is (i) For a given axial ratio, and 
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perfect alignment, find the aligned grain size distribution by fitting the frequency 
dependence of interstellar polarization. Compare this to a consistent model of the 
frequency dependence of interstellar extinction, keeping track of the mass of all 
components (unaligned grains contribute to t and not p, and so cause dilution). 
For this we use models by Kim & Martin l)1995|l : only the silicates are aspherical 
and aligned, which is certainly adequate, (ii) Calculate {p/T)ex for this model (it 
can be evaluated for frequencies throughout the optical). We took the average 
of the two estimates in Table 1 of Kim & Martin I)1HH5(I and noted the implied 
substantial systematic error, (iii) Compare this to the observed maximum (p/r)ea; 
(Fig.[7I) to deduce a reduction factor R^x < 1 reflecting the imperfect ahgnment. 
(iv) For the same model (shape, axial ratio, grain components), self-consistently 
calculate the polarization of the low-frequency thermal emission of the aligned 
grains, pUem (e.g., Fig^. (v) Assuming that Rem = Rex^ apply R from the in- 
terstellar polarization model for that axial ratio to find the expected maximum 
degree of polarization Rpa^m for the disaligned grains. 

How robust this estimate is can be judged by repeating the recipe calculation 
for different shapes and axial ratios, (vi) There is a further complication, dilu- 
tion of the polarization because of thermal radiation by non-aligned grains. This 
reduction factor d < 1 needs to be evaluated as self-consistently as possible (§|OJ. 

5 Results and Discussion 

The results from the different steps of this recipe are collected in Tabic [21 There 
puem is evaluated at 350 GHz, but this choice is inconsequential (Fig. EJ. 

For perfectly aligned oblate silicate particles (i? ~ 1 in this case), the axial ratio 
needs to be no higher than 1.4 to produce the maximum {p/t)ex observed. For 



Table 2. Prediction of the Net Submillimctre Polarization 



Oblate Prolate 



Axial ratio 


V2 


2 


4 


6 


21 


22 


42 


{P/T)ex 


0.04 


0.08 


0.12 


0.12 


0.05 


0.09 


0.14 


R 


0.67 


0.36 


0.23 


0.23 


0.53 


0.30 


0.19 


P^em 


0.31 


0.56 


0.83 


0.90 


0.39 


0.56 


0.91 


Rpaem 


0.21 


0.20 


0.19 


0.21 


0.21 


0.17 


0.17 


dRpaem (%) 


9.3 


8.9 


8.7 


9.4 


9.3 


7.5 


7.8 



1 perfect spinning alignment 

2 picket fence alignment 
^ evaluated at 350 GHz 
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larger axial ratios, grains must be somewhat disaligned by a quantifiable amount 
(i? < 1) to produce the same {p/T)f.x- The value oipa^m is large, and as expected 
higher for the more extreme axial ratios (Fig. ^J. The application of the recipe 
to the models for axial ratios closer to unity is probably more reliable, and these 
ratios are possibly the most realistic as well. But it turns out that the product 
Rpaem is fairly robust, depending little on the shape and axial ratio. 

5.1 Dilution 

Non-aligned grains contribute to the thermal emission and dilute the maximum 
polarization expected. The Kim & Martin models used involve silicates and (large) 
graphite particles (Kim et al. I1994|l . following Mathis et al. H1977|l and Draine 
& Lee (1984). The relative contributions to the submillimetre emission can be 
judged from calculations by Draine & Anderson H1985|l . which do produce about 
the right amount of far-infrared emission per unit column density. This suggests 
quite a reduction, by o? ^ 1/3. A more recent variant of this model, with slightly 
different size distributions and apportionment of material (Li and Draine 1200 l|l . 
gives d ~ 0.58 at 353 GHz (their Fig. 9). We adopt the mean, d ~ 0.45, and note 
another major contribution to the systematic uncertainty. 

Note that the dilution of {p/T)ex in the optical by extinction by unaligned 
graphite (carbonaceous) particles has resulted in a larger Rex, and so this low 
d can be considered as payback, compensating in the product dR. Likewise, if 
there were optical polarization by large carbonaceous particles, there would be an 
accompanying decrease in R but an increase in d, and so this change too would 
not produce a greatly different net Pem in the submillimetre. 

5.2 Predicted Net Polarization pem 

This self-consistent model then predicts a maximum net polarization dRpUem as 
listed in Table|21and summarized as the median pem(%) = 8.9±0.7±3.5. The result 
is not very dependent on the shape or axial ratio of the grains, as reflected in the 
modest rms 0.7% among the different models. The systematic uncertainty 3.5% 
is much more substantial, with the rough estimate based on actual application 
of the recipe to the models: calculating Tex appropriately (step (ii); S 14. 1() . the 
(frequency dependent) dilution d near 353 GHz (step (vi); S 15. 1(1 . and encompassing 
the assumption Rem = Rex (step (v); S I4.1|I . 

Thus the maximum polarization to be expected is quite large, even without the 
possible boost accompanying spectral flattening of the amorphous silicate emissiv- 
ity through 350 GHz. 

But when averaged over large regions with (i) non-uniform alignment (beam 
dilution) or (ii) less than the optimal alignment, including the effects changes in 
the direction of the magnetic fleld, or unfavorable field orientation (something that 
changes systematically on a large scale in the Galaxy), or (iii) alignment which 
changes along the line of sight, then typically half of this might be expected, 
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judging from Fig.|Slwhich shows the depolarization from (ii) and (iii). This would 
still leave p^m ^5%. 

5.3 Observations 

Benoit et al. ()2U(J4|I present Archeops measurements of polarization in the Galactic 
plane at 353 GHz. They find Pem ~ 4 — 5% for the diffuse emission averaged over 
several square degrees and even higher values in some large clouds. The orientation 
of the E vector roughly perpendicular to the plane is reassuringly as expected. 
However, no clear correlation with optical interstellar polarization (for which the 
path lengths to the background stars would be shorter than for the emission) has 
been established and the data are not of high enough signal to noise to follow this 
at high resolution. Because of the long lines of sight and high column densities, 
the Galactic plane, even its diffuse emission, is probably not the best place to 
evaluate the polarization of the higher latitude cirrus foreground of the CMB. 
The amount of polarization observed is actually close to what is predicted by the 
analysis above. Although such estimates were actually made in 2003, before the 
Archeops results, this can hardly be taken as a triumph, given the large systematic 
uncertainties and the extra uncertainties of modeling the emission in the Galactic 
plane. 

Nevertheless it does support the view that the polarization of the cirrus at 
higher latitude will be significant. Although the variations in alignment along and 
across the line of sight at high latitude are not known, it seems possible (because of 
the shorter paths and thus simpler geometry being integrated) that these variations 
will be less than in the Galactic plane, resulting in less depolarization. 

5.4 Spectral Behaviour ofpem 3nd P 

As noted above, the spectral index of the amorphous silicate emissivity might be 
less steep than the f3 ^ 2 inherent in our calculations. A modest change was 
introduced by Li & Draine |2OT1] their eq. 1 and Fig. 9) to fit the frequency 
dependence of the diffuse emission. Boudet et al. 1)2005(1 suggest even stronger 
variations with frequency over the range 150 - 3000 GHz. In this case, the factor 
d would be tend to increase with decreasing frequency as the relative importance of 
large silicates to the submillimetre emission rose. In models where it is the silicates 
that produce the polarization (subscribed to here) , this frequency dependence of d 
would result in an increase in the net Pem at the lower frequencies and our estimate 
at 353 GHz could be low by a factor 1.6 (not included in the systematic errors 
below). It will be interesting to learn the results from B2K (and Planck) which 
measure within this interesting range of frequencies. 

Because of different weighting of different components, the spectral dependence 
of the polarized intensity P can be different than that of total emission, / (§12). 
Through dilution, non-aligned grains lower the net polarization. If the frequency 
dependence of emission for the diluting component is different, then this introduces 
a frequency dependence into Pem = P/P We have just mentioned one example. 
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To characterize the intrinsic emission properties of the carriers of the polar- 
ization more directly, this confusion from dilution can be avoided in principle by 
examining P, which would isolate the spectral dependence of the polarizing emit- 
ters. For frequency ranges over which paem for this component could be assumed 
to be constant, this would directly give (3, and its variations, for the spectral 
emissivity (/) of this aligned-grain component. 

P is also the quantity on which to base statistical evaluation of the spatial 
variations. 

A toy model for illustration is presented in FigureEl Here it has been assumed 
that (3 for the aligned silicate component changes significantly with frequency, in 
the spirit of what has been inferred from astronomical observations and measur- 
ments in the laboratory (Boudet et al. 120051 J. -P. Bernard, private communica- 
tion) , although to be sure this is not (yet) well constrained. Specifically the model 
uses P = 1.7 for A < 1 mm {v > 300 GHz), flattening (perhaps too abruptly in 
this illustrative model) to 0.5 through 2 mm (150 GHz), and then steepening to 
1.5 (again for visual emphasis) beyond 5 mm {ly < 60 GHz). For the silicate grains 
we take T = 17 K and puem = 0.2 (0.089/0.45 as in 15.21 In a more realistic 
model there would be a slight increase in paem with decreasing frequency as a 
resultof impliedchanges in the complex refractive index of the silicate. There is 
also an unpolarized graphite component with T = 20.5 K with relative emission 
such that d — 0.45 at 353 GHz (consistent with Table[21). For the graphite compo- 




100 1000 

(GHz) 



Fig. 9. Frequency dependence of observables pem, I, and P (solid curves, from top) in 
an illustrative model in which the value of /3 for the amorphous silicate component varies 
with frequency (see text) . The silicate and graphite constributions to I are shown by the 
dashed and dot-dash curves, respectively (for normalization, see text). 
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nent, (3 = 2. It should be kept in mind that the "graphite" might alternatively be 
amorphous carbon, in which some frequency dependent changes in [3 might also 
be expected. 

The "observables" , /, P and pem, are shown by the sohd curves, in FigurelHl 
The model illustrates several important features. The spectral index measured 
from J is a function of frequency and would not be representative of either grain 
material. On the other hand, the spectral index measured from P would track 
that of the aligned amorphous silicates, even though this component is diluted in 
/. Toward lower frequencies silicates provide a more dominant contribution to the 
total emission {d is larger) and so even though paem does not rise significantly for 
this component, the net Pem does increase to lower frequencies. 

It appears that the frequency dependence could be quite complicated, such that 
it would not be correct to extrapolate at constant (3 and/or p^m from meaurements 
at a single high frequency toward the CMB frequncies. Rather, to disentangle this 
most convincingly, multi-frequency measurements of both pem and P (and /) will 
need to be examined in regions of relatively bright foregrounds to produce a self- 
consistent model. Such an assessment of the dust polarization on the brighter 
cirrus, uncontaminated by the CMB at lower frequencies, would provide a basis 
for understanding the fainter cirrus contaminating the CMB. While this is not 
possible with Archeops, B2K and Planck have polarization capability at three 
frequencies and PILOT will extend this to higher frequencies (Table^). It would be 
advantageous to make use of polarization measurements below 100 GHz too (e.g., 
WMAP and Planck), but at some point synchrotron emission, which probably has 
a higher intrinsic polarization, would no longer be negligible (Fig.^, and there is 
the anomalous emission to considere too.. 

5.5 Alignment of Small Grains and Polarization of Anomalous Emission 

As discussed in S 13.21 {p/t)qx is very low in the ultraviolet, where the extinction 
comes from small grains (VSGs, PAHs). What polarization there is (Fig. O is 
consistent with a fading contribution coming from big grains (Kim & Martin [T995|l . 
However, it is possible that there is a very low level of residual alignment of small 
grains, such as from the Davis-Greenstein process, could be present (Kim & Martin 
[BMI Wolff et al. nWn Lazarian & Finkbeiner IM)^ . 

Tiny grains, including if not exclusively PAHs, are the ones that could spin 
most rapidly to produce low frequency (~ 20 GHz) emission (Draine & Lazarian 
I1998|l . possibly accounting for the dust-correlated anomalous microwave emission 
(Finkbeiner "^IM", Finkbeiner et al. OIMl de Oliveria-Costa et al. IMTl Davies 
12006 . Davis 2006). Since these small particles should not be well aligned, the 
degree of polarization of that component would not be high (tiny compared to the 
degree of polarization of the synchrotron component), a diagnostic feature that 
might be used to advantage by WMAP and Planck LFI. Still, the spectre of a few 
percent polarization for spinning dust could be a serious contamination for precise 
CMB polarization measurments (Lazarian & Draine I^OOOI Lazarian & Finkbeiner 
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We can be somewhat more quantitative in a semi-empirical if model-dependent 
way, starting from the Li & Draine I)2U(J1|I model, in which PAHs are responsible 
for most if of the strong 2175 A feature in the extinction curve (as reasonable and 
economical assumption) . We then use the tight limits placed on the polarization of 
the 2175 A feature (Martin et aL 119991 ^ I3.2|) . The relative change in polarization, 
compared to extinction, at the bump gives for this component p/t^x < 0.002. The 
same quantity for perfectly aligned PAHs could be calculated using small-particle 
formulae as in equations 12.31 and 12.41 (the extinction is pure absorption) , given a 
shape and refractive indices. However, there are no refractive indices since only 
absorptivity has been modeled (Li & Draine l?001|l . and for these large molecules 
it might not even be particularly appropriate. Nevertheless, since the PAHs are 
probably highly flattened, this quantity would be near unity, just as it is for tiny 
graphite spheroids of even modest asphericity (Martin et al. 119951 Wolff et al. 
I1997|l . Therefore, Rex < 0.002. Similarly, for the rotational dipole emission of 
well aligned tiny grains, paem would also be unity. It would be reasonable to 
assume similar effects from disorientation: Rem ~ Rex- Thus for the dust-related 
anomalous emission, we would predict a degree of polarization less than 0.2%, 
probably very much less given that polarization of the 2175 A feature is unusual. 
This is more than an order of magnitude, probably nearer two, less than for the 
degree of polarization of the thermal dust emission, and also very much less than 
expected for synchrotron emission. The real situation is more subtle, however, 
because polarization of the rotational emission depends on the alignment of the 
angular momentum vector with respect to the magnetic field rather than the body 
axis with respect to the field, needed for the ultraviolet polarization (Lazarian & 
Draine IMTHjl . If there is only partial alignment of the angular momentum vector 
and the body axis, then the low frequency polarization would be higher and with 
detection of polarization of the anomalous emission one might learn something 
more about the alignment. 

To complicate things further, there is an alternative model for the anomalous 
microwave emission, magnetic-dipole emission from magnetic grains (Draine & 
Lazarian 1 19991 Lazarian & Finkbeiner I2003|l . which would have a very distinctive 
frequency-dependent polarization signature. Not enough is known about the prop- 
erties of these grains or their alignment to make firm predictions of the degree of 
polarization. However, since magnetic inclusions might be a factor in grain align- 
ment, this emission might be more closely related to "normal-sized" grains which 
are better aligned, and therefore might be more highly polarized than emission 
from spinning dust. Following similar bootstrapping arguments such as above, the 
order of magnitude might be like that given by Rpaem in Table |2 thus ~ 20%, 
reduced by other diluting contributions to the microwave emission. Again, syn- 
chrotron polarization would be a consideration. 
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Appendix: The Last Word — Grain Alignment Explained 

I recommend looking at the above-mentioned Planck web pages, but one tiling that cannot be 
found there is an explanation of why the grains are aligned. Recall that somehow the dust 

particles are supposed to be spinning preferentially around the magnetic field. At the conference 
I announced my forthcoming monograph "Alignment for Dummies" (coming soon to a discerning 
supermarket checkout counter near you) and to pique curiosity further (advanced sales benefit 
from feverish anticipation) 1 elaborated on the merits of one new theory of grain alignment. It 
turns out that the new theory is rather straightforward, so much so that 1 attempted to explain 
it en frangais, hoping to motivate, indeed enable, Lagache and colleagues to update their Planck 
web pages. 

C'est facile, comme "un, deux, trois. " 

Un! Voila, une grande poussiere. Je pense que c'est organique, d la Greenberg. 
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Fig. 10. A moment of great anticipation in the demonstration of "Intelligent Alignment," 
at the point of "contact" as the propeller motion is energized. 

Deux! Voila, le champs magnetique interstellaire - plus grand que les Champs- Ely see. 
Mais, aucun alignement! Dommage! 
Trois! Voila, la main de Dieu! 

Alors, attention! Un: une poussiere. Deux: un champs magnetique. Trois: la main de Dieu. 
Alignement! 

The live demonstration (Fig. llOi does benefit from a grocery bag containing appropriate 
ingredients! Any good theory needs a name and it is suggested that this theory might well merit 
the status of "Intelligent Alignment." 



